Hidden dimer order in the quantum compass model 
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We introduce an exact spin transformation that maps frustrated ZijZij+i and XijXi+ij spin 
interactions along the rows and columns of the quantum compass model (QCM) on an L x L square 
lattice to (L — 1) x (L — 1) quantum spin models with 2(L — 1) classical spins. Using the symme- 
try properties we unravel the hidden dimer order in the QCM, with equal two-dimer correlations 
(XijXi+ijXkjXk+xj) and (Xi^Xi+i^Xi^Xi+i^) in the ground state, which is independent of the 
actual interactions. This order coexists with Ising-like spin correlations which decay with distance. 
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The quantum compass model (QCM) originates from 
the frustrated (Kugel-Khomskii) superexchange^ in tran- 
sition metal oxides with degenerate 3d orbitals. Recent 
interest in this model is motivated by its interdisciplinary 
character as it plays a role in the variety of phenomena 
beyond the correlated oxides. It describes a quantum 
phase transition between competing types of order when 
anisotropic interactions are varied through the isotropic 
point, as shown by an analytical method^ mean- field 
(MF)^ and numerical studies^ The QCM is dual to 
the models of p + ip superconducting arrays^ and to 
the toric code model in a transverse field It was also 
suggested as an effective model for Josephson arrays of 
protected qubits^ as realized in recent experiment^ Fi- 
nally, it could describe polar molecules in optical lattices 
and systems of trapped ionsJ^ 

In spite of several numerical studies^ - — the nature of 
spin order in the two-dimensional (2D) QCM is not yet 
fully understood. By an exact solution of the QCM on a 
ladder we have shown, however, that the invariant sub- 
spaces may be deduced using the symmetry^ The 2D 
QCM shows a self-duality^ which might serve to reveal 
nontrivial hidden symmetries^ In this Letter we employ 
exact spin transformations which allow us to discover a 
surprising hidden dimer order in the QCM which mani- 
fest itself by exact relations between four-point correla- 
tion functions in the ground state. We also demonstrate 
nonlocal MF splitting of the QCM in the ground sub- 
space and determine spatial decay of spin correlations in 
the thermodynamic limit. 

Reduced Hamiltonian. — We consider the anisotropic 
ferromagnetic QCM for pseudospins 1/2 on an L X L 
square lattice with periodic boundary conditions (PBC), 



i,j=l 

where {X^j,Z^j} stand for Pauli matrices at site 



i.e., Xij = af j and Zij = o~~j components, interacting 
on vertical and horizontal bonds. In case of L being even, 
this model is equivalent to the antifcrromagnctic QCM. 



We can easily construct a set of 2L operators which com- 
mute with the Hamiltonian but anticommutc with one 



another £ p 



nLi 



Xij and Qj = Yl i= i Zij- Below 
we will use as symmetry operations all Ri = RiPi + \ and 
Qj to reduce the Hilbert space; this approach led to the 
exact solution of the compass ladder^ The QCM Eq. 
(1) can be written in common eigenbasis of {Ri, Qj} op- 
erators using: 
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where Zqj = 1 and X' iQ = 1. One finds that the trans- 
formed Hamiltonian, T-L'(a) = —(1 — ct)H' x — aH' z , con- 
tains no Xj,.j and no Z[ L operators so the corresponding 
Zj^ j and X[ L can be replaced by their eigenvalues qj and 
ri, respectively. The Hamiltonian ri'(a) is dual to the 
QCM H(a) of Eq. (1) in the thermodynamic limit; we 
give here an explicit form of its x-part, 
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where r = Yli=i r i> an< ^ new 

nonlocal Pj = JJ^J X' p/J 
operators originate from the PBC. The z-part H' z follows 

Z[ p and by 



from H' x by lattice transposition X- j 
ri — > Sj = qjqj+\. Ising variables and Sj are the eigen- 
values of the symmetry operators Ri and Sj = QjQj+\. 

Instead of the initial L x L lattice of quantum spins, 
one finds here (L — 1) x (L — 1) internal quantum spins 
with 2(L — 1) classical boundary spins. The missing spin 
is related to the Z2 symmetry of the QCM and makes 
every energy level at least doubly degenerate. Although 
the form of Eq. (0| is complex, the size of the Hilbert 
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space is reduced in a dramatic way by a factor 2 2i_1 — 
which makes it possible to perform easily exact (Lanczos) 
diagonalization of 2D Lx L clusters up to L = 5 (L = 6). 

Equivalent subspaces. — The original QCM of Eq. (1) 
is invariant under the transformation X' f-> Z', if one 
also transforms the interactions, a O (1 — a). This im- 
plies that subspaces (r, s) and (s, r) give the same energy 
spectrum which sets an equivalence relation between the 
subspaces — two subspaces are equivalent means that the 
QCM (1) has in them the same energy spectrum. This 
relation becomes especially simple for a = \ when for all 
ri's and Si's subspaces (r, s) and (s, r) are equivalent. 

Now we will explore another important symmetry of 
the 2D compass model reducing the number of nonequiv- 
alent subspaces — the translational symmetry. We note 
from Eq. (TjTJ) that the reduced Hamiltonians are not 
translationally invariant for any choice of (r , s) even 
though the original Hamiltonian is. This means that 
translational symmetry must impose some equivalence 
conditions among subspace labels {r, s}. To derive them, 
let's focus on translation along the rows of the lat- 
tice by one lattice constant. Such translation does not 
affect the Pi symmetry operators, because they con- 
sist of spin operators multiplied along the rows, but 
changes Qj into Qj+i for all j < L and Ql — > Qi- 
This implies that two subspaces (f, qi, q^, . . . , qi,) and 
(r, ql, qi, qi, . . . , qi,-i) are equivalent for all values of r 
and q. Now this result must be translated into the lan- 
guage of (r,s) labels, with Sj = qjqj+i for all j < L. 
This is two-to-one mapping because for any s one has 
two <fs such that q + = (1, s\, S1S2, • • • , S1S2 . . . Si-i) 
and g_ = —q+ differ by global inversion. This 
sets additional equivalence condition for subspace la- 
bels (r,s): two subspaces (r, it) and (r,v) are equiv- 
alent if two strings (1, u\, U1U2, . . . , U1U2 ■ ■ ■ Wl— 1) and 
(1, vi,v\V2, ■ ■ ■ , i>i«2 ■ ■ ■ vl-i) are related by translations 
or by a global inversion. For convenience let us call these 
two vectors TI (translation-inversion) related. Lattice 
translations along the columns set the same equivalence 
condition for r labels. Thus full equivalence conditions 
for subspace labels of the QCM are: 

• For a = I two subspaces (r, s) and (it, v) are equiv- 
alent if r is TFrelated with u and s with v or if r 
is Tl-rclated with v and s with u. 

• For a ^ I two subspaces (r, s) and (it, v) are equiv- 
alent if r is Tl-related with u and s with v. 

We have verified that no other equivalence conditions ex- 
ist between the subspaces by numerical Lanczos diagonal- 
izations for lattices of sizes up to 6 x 6, so we can change 
all if statements above into if and only if ones. 

Hidden dimer order. — Due to the symmetries of the 
QCM Eq. (1) only {Z hJ Zij +d } and (X itj X i+ d,j) spin 
correlations are finite (d > 0). This suggests that the en- 
tire spin order concerns pairs of spins from one row (col- 
umn) which could be characterized by four-point corre- 
lation functions of the dimcr-dimer type. Indeed, exam- 
ining such quantities for finite QCM clusters via Lanczos 
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FIG. 1: (color online). Example of the application of the 
proved identities in two cases: left panel — Eq. (|10[) for a 
chosen dimer Z\j = XijXi+ij (label 1), correlations of T>i.j 
with a given Z\,; are the same for dimers T>k,i marked with 
the same label (color); right panel — Eq. (jlll) long range 
correlation function (XijXi+j+ij) along the column (circles) 
is equal to the 2d-point (XX . . . X) correlation function along 
the row (frame of length d) . 



diagonalization we observed certain surprising symme- 
try: for any a dimer-dimer correlators (DijDk i), with 
Vij = Xi jXi+ij, are invariant under the reflection of 
the second dimer with respect to the diagonal passing 
through site (i,j), see left panel of Fig. [T] This general 
relation between correlation functions of the QCM will 
be proved below. 

We will prove that in the ground state of the QCM for 
any two sites and (fc, I) and for any < a < 1: 

(Xi,jXi+ijXk,iXk+i,i) 

= (XijXi + ijXi-5 t k+8Xi-S+l,k+8), (5) 

where 5 = j — i,™ i.e., the second dimer is reflected with 
respect to the diagonal. To prove Eq. §5§ let us trans- 
form again the effective Hamiltonian Q in the ground 
subspace (r, = Sj = 1) introducing new spin operators 

j 

= Zi,jZi t j + i, X-j = Y[ Xi,r, (6) 

r=l 

with i,j = 1, . . . , L — 1 and Zi l = L This yields to 

i-lL-l L-li-l 
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I Z a ^ + (^Z a ,iZ a .i + i + Zi ta Zi±x ia J > 

a { b i=l ) 



L-2 L-2 

E E ZijZi,j+iZi+i,jZi+i,j+i, 
1=1 j=i 
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where a, b = 1,L — 1. Due to the spin transformations 
(|2I3I6I) . Xij operators are related to the original bond 
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FIG. 2: (color online). Schematic view of the z-part of the 
reduced ground subspace Hamiltonian (j8}: circles in the cor- 
ners stand for Zi,j spin operators related to the site (i, j), 
dashed (red) frames are ZZ operator products acting on the 
boundaries of the lattice, and solid (blue) square stands for 
one of the plaquette ZZZZ spin operators. The exemplary 
three sites in the identity (|10[1 are: (i, i), (k,l) and (l,k). 



operators by XijX i+ ij — X,ij, which implies that 



(9) 



Because of the PBC, all original Xij spins are equiv- 
alent, so we choose i = j. The x-part ([7|) of the Hamil- 
tonian is completely isotropic. Note that the z-part © 
would also be isotropic without the boundary terms (see 
Fig. [2]); the effective Hamiltonian in the ground sub- 
space has the symmetry of a square. Knowing that in 
the ground state we have only Zi degeneracy, one finds 



(Xi^Xkj) = (X iyi Xi t k), 



(10) 



for any i and (k, I). This proves the identity ([5]) for S = 0; 
(5/0 case follows from lattice translations along rows. 

The nontrivial consequences of Eq. (fT0| are: (i) hidden 
dimer order in the ground state of the QCM, i.e., an 
" isotropic" behavior of the two-pair correlator in spite of 
anisotropic interactions in the entire range of < a < 1 
(see Fig. [3]), and (ii) long range two-site (XijX i+c t+i,j) 
correlations of range d along the columns which are equal 
to the multi-site (XX ... X) correlations involving two 
neighboring rows, see right panel of Fig. [1] The latter 
follows from the symmetry properties of the transformed 
Hamiltonian (|7|8|) applied to the multi-site correlations: 
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Mean-field approximation. — The x-part of the Hamil- 
tonian obtained from Eq. (£[]) in case of open boundaries 
reads: 
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jX'i.j+i + X[ x + TiX[ 
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FIG. 3: (color online). Two-dimer (Xl, 1X2,1X^,1X^+1,1) cor- 
relations for L = 6 and 0.2 < a < 0.8: (k, I) = (1,2), (1,3) 
and (1,4) are shown by solid, dashed and dotted line. 



and similarly for the z-part. In the ground subspace 
(rj = 1) this resembles the original QCM Eq. (1) but 
with linear boundary terms, which should not affect the 
ground state properties in the thermodynamic limit and 
can be regarded as symmetry breaking fields, resulting in 
finite values of (X- j ) and (Z[ ^ ) . Omitting the boundary 
terms in H' x and H' z and putting infinite L we recover 
the 2D QCM written in nonlocal primed spin operators. 
Now we can construct a MF splitting of the 2D lattice 
into (ferromagnetic) Ising chains in transverse field, tak- 
ing (Z 1 ) = (Z'ij) as a Weiss field for each row i: 



E 
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{ (1 - a)Xl 3 Xl 3+l + 2a(Z')Z' id } . (13) 



In analogy to the compass ladder^ it can be solved by 
Jordan- Wigncr transformation for each i: 



X 'iJ 



1 



2cL .< 



i,3 »,J ' 



c- e 4 



(14) 

U(l-2cl r c iir ), (15) 



introducing fermion operators {cjj}. The diagonaliza- 
tion of the free fermion Hamiltonian can be completed 
by performing first a Fourier transformation (from {j} to 
{k}) and next a Bogoliubov transformation (for k > 0): 

7fe = «fc4 + Pt c -k and 1-k = a k c l + Pk c -k' whcrG 
{&t,l3t} are eigenmodes of the Bogoliubov-de Gennes 
equation for the eigenvalues ±Ek (with Ek > 0). The re- 
sulting ground state is a vacuum of 7! fermion operators: 

*o) = Ilfc>o( a fc + ^ c -fe c l)l )' wmcl1 can serve to cal- 
culate correlations and the order parameter of the QCM 
in the MF approach. In agreement with numerical results 
(not shown), the only nonzero long range two-site spin 
correlation functions are: (XijXi+d.j) and (ZijZi_j+d). 
For d > 1 they can be represented as follows: 



(X i}j X i+d ^) - {X- .xl , 



(ZijZij+d) 



( Z i,j Z i.3 + l 



Z i,j+d- 



(16) 
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FIG. 4: (color online). Long range spin correlations of the 2D 
QCM Eq. (1) obtained in the MF approach for a > ~. Lines 
starting from 1 at a = 1 are the (Zi^Zij+j) correlations (|17p 
for d = 1,2,3,4,5,10,20,80, while lines starting from at 
a = 1 are the (XijXi+dj) correlations (|16[1 with d = 1,2,3; 
both in descending order. 

Having solved the self-consistency equation for (Z') = 
(1 - 2<n», with (n) = iE M K 2 + 4 +2 ). one can 
easily obtain (X- jX^j +1 ) (JTHJ) for increasing a: 

fe>0 

+ sinfc(a^- -«+/?+)}. (18) 

The nonlocal ^Z' . . . correlations (jTTJ) are 

more difficult to find but they can be approximated by 

(Z[,Z' itj+1 . . . Z[, +d _ x ) = n |«+ 2 (l - 2 £) % /3+ 2 j , 

(19) 

where i -> oo and k = (21 - l)f with I = 1, 2, . . . , -§. 
This approximation is valid as long as rf < L. One finds 
that the long range (ZijZij +c i) correlations in ^-ordered 
phase at cv > i show the absence of the Ising-like long 
range order for a < 1 (Fig. HJ) — they decrease slowly 
with growing distance d or decreasing In contrast, 



the (Xi,jXi + d,j) correlations arc significant only for near- 
est neighbors (d = 1) and close to a = i. 

The advantage of this nonlocal MF approach for the 
QCM Eq. (1) over the standard one, which takes (Z) 
as a Weiss field, is that we do not break the {Pi,Qj} 
and Z2 symmetries of the model. What more, thanks 
to numerical and analytical results we know that order 
parameter of the QCM is given by (H z )^ — the quantity 
behaving more like (Z') rather than (Z) (having (Z) > 
would mean long range magnetic order !). Another inter- 
esting feature of the Hamiltonian (fTJ) is that it describes 
all nonlocal compass excitations over the ground state, 
while the local ones manifest themselves by directions of 
symmetry breaking fields. These nonlocal column (row) 
flips are especially interesting from the point of view of 
topological quantum computing 2 - because they guarantee 
that the system is protected against local perturbations. 

Conclusions. — On the example of the QCM, we argue 
that the properties of spin models which are not SU(2) 
symmetric can be uniquely determined by discrete sym- 
metries like parity. In this case conservation of spin par- 
ities in rows and columns, for x and z-components of 
spins, makes the system in the ground state behave ac- 
cording to a nonlocal Hamiltonian In the ground 
state most of the two-site spin correlations vanish and 
the two-dimer correlations exhibit the nontrivial hidden 
order. The excitations involve whole lines of spins in 
the lattice and occur in invariant subspaces which can 
be classified by lattice translations — the reduction of 
the Hilbcrt space achieved in this way is important for 
future numerical studies of the QCM and will play a 
role for spin models with similar symmetries. Finally, 
the nonlocal Hamiltonian containing symmetry break- 
ing terms suggests the MF splitting respecting conser- 
vation of parity and leading to the known physics of one- 
dimensional quantum Ising model describing correlation 
functions and the order parameter of the QCM. 
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The symmetry of Eq. (1) implies a similar relation for Z- 

correlations with dimers along the horizontal bonds. 

At a = \ the Z-ordered and X-ordered Ising-like phases 



of the QCM are degenerate, see e.g. Ref. @. 

The nonlocal Hamiltonian is just the original QCM Eq. (1) 

written in eigenbasis of {Pi,Qj} parities. 



